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1.  INTRODUCTION 


The  following  is  an  extension  of  an  earlier  theoretical  study  of  the  H  +  OCS  potential  energy  surface 
(PES)  (Rice,  CarUand,  and  Chabalowski  1993).  In  that  study,  qualitative  features  of  the  following  two 
reaction  channels  were  calculated  at  the  PUMP4/UMP2/6-31G'''*  level. 


H(2s)  +  OCS(^I)  ->  OHf^n)  +  CS(^I) 

AH°298  "  kcal/mol 

(I) 

H(2s)  -t-  OCS('S)  SfR^n)  +  CO(^I) 

AH®298  “  “12.1  kcal/mol 

(11) 

In  this  study,  we  attempt  a  more  quantitative  prediction  for  the  critical  points  on  this  surface.  The 
calculations  presented  here  use  basis  sets  as  large  as  Dunning’s  aug-cc-pVTZ  (Dunning  1989;  Kendall, 
Dunning,  and  Harrison  1992;  Woon  and  Dunning  1993)  and  correlation  corrections  as  high  as  PUMP4 
and  QCISD(T)  (Poplc,  Hcad-Gordon.  and  Raghavachari  1987). 

In  the  previous  study  (Rice,  Cartland,  and  Chabalowski  1993),  we  found  structures  corresponding  to 
minima  and  saddle  points  along  die  reaction  pathways  for  (I)  and  (II).  These  consisted  of  6  stable  4-body 
complexes  and  12  saddle  points,  enabling  us  to  characterize  the  primary  reaction  paths  for  (I)  and  (II). 
The  study  showed  that  the  low-energy  approach  of  the  hydrogen  to  OCS  is  broadside,  rather  than  end-on 
for  both  reactions.  The  study  substantiated  experimental  hypotheses  of  the  existence  of  stable  four-body 
reaction  intermediates  (BOhmer,  Mikhaylichenko,  and  Wiltig  1993;  Nickoiai.scn  and  Cartland  1993; 
Hausler,  Rice,  and  Wittig  1987),  and  the  "tight"  four-body  transition  slates  leading  to  the  products  of  (II) 
(Tsunashima  ct  al.  1975;  Lee,  Stief,  and  Timmons  1977).  Additionally,  we  offered  qualitative  explanations 
for  the  observed  noastatistical  product  energies  di.stributions  for  Reaction  (II)  and  the  statistical  product 
energy  distributions  for  Reaction  (I)  (BOhmer,  Mikhaylichenko,  and  Wittig  1993;  Nickolaisen  and  Cartland 
1993;  Hausler,  Rice,  and  Wittig  1987).  The  results  of  the  work  presented  here  do  not  change  those 
explanations  or  conclusions.  Unfortunately,  the  level  of  theory  and  basis  set  used  in  oui  previous  study 
prevented  us  from  making  quantitative  comparison  of  energetics  with  measured  activation  energies  for 
Reaction  II  (I’sunashima  ct  al.  1975;  Lee,  Stief,  and  Timmons  1977). 

The  thermal  rate  measurements  of  Tsunashima  et  al.  (1975)  and  Lee,  Stief,  and  Timmons  (1977) 
provided  Arrhenius  parameters  for  Reaction  (II).  The  activation  energy  and  pre-exponential  factor  reported 
by  Tsunashima  ct  al.  (1975),  for  rates  measured  between  300  to  525  K,  are  3.90  ±0.370  kcal/mol  and 
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1.5  ±0.20  X  10“''  cm^^molcciile-s.  Lee,  Stief,  and  Timmons  (1977)  report  similar  Arrhenius  parameters 
for  measurements  between  261  to  500  K;  they  give  an  activation  energy  and  pre-exponential  factor  of 
3.85  ±0.110  kcal/mol  and  9.06  ±1.53  x  10“'^  cmVmolccule-s,  respectively.  The  magnitudes  of  the 
reported  pre-exponential  factors  are  small  compared  to  other  hydrogen  abstraction  reactions,  leading  die 
experimentalists  to  conclude  that  u  light-activated  complex  was  involved  in  the  formation  of  product 
(Tsunashima  et  al.  1975;  Lee,  Stief,  and  Timmons  1977).  Our  previous  results  supported  this  conclusion 
of  a  tighi-activalcd  complex.  Our  calculations,  however,  had  enough  error  in  the  energetics  dial  the 
activation  energy  issue  could  not  be  addressed. 

We  present  results  of  extended  basis  set  calculations  at  various  levels  of  correlation  correction  for 
stationary  points  on  the  H  +  OCS  PES.  The  purpo.se  of  this  study  is  to  provide  a  quantitative  comparison 
of  these  energies  to  the  experimental  results,  and  to  substantiate  the  existence  and  structure  of  the  critical 
points  located  in  the  earlier  study  (Rice,  Caitland,  and  Chabalowski  1993).  We  have  also  attempted  to 
assess  the  accuracy  of  these  calculations  in  die  comparison  with  experimental  data.  As  die  quality  of  the 
theorcUcal  treatment  improves,  a  trend  in  the  predicted  energies  emerges  that  shows  convergence  toward 
the  experimental  values.  This  allows  for  a  rough  estimate  of  the  uncertainty  in  die  energies  calculated  for 
the  structures  on  the  PES. 

A  question  wc  did  not  address  in  the  previous  study  was  the  role  of  hydrogen  migration  in  the 
reactions.  Besides  locating  each  stationary  point  that  was  found  in  our  previous  study  at  the  lower  level 
of  theory  and  basis  set,  wc  have  located  the  transidon  state  for  the  hydrogen-migration  reaction.  The  size 
of  this  barrier,  however,  makes  this  a  high-energy  reaction. 

2.  METHODS 

The  calculations  presented  here  were  performed  using  the  Gaussian  92  set  of  codes  (Frisch  ct  al. 
1992).  All  geometry  optimizations  meet  or  exceed  the  default  values  for  convergence  in  diesc  codes.  No 
attempts  were  made  to  correct  for  basis  set  superposition  error  (BSSE)  (Boys  and  Bcmardi  1970; 
Andzclm,  Klobukowski,  and  Radzio-Andzchn  1984)  since  studies  at  the  SCF  level  on  the  interaction 
energy  of  HE  ±  HE  for  34  different  basis  sets  (Schwenkc  and  Truhlar  1985,  1986)  showed  that  the 
counterpoise  method  (Boys  and  Bcmardi  1970),  or  variations  of  this  method,  gave  a  corrected  energy  tliat 
"was  no  more  reliably  accurate  than  the  uncorrcclcd  energy"  (Davidson  and  Feller  1986).  MP  calculations 
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reported  here  were  run  without  any  frozen  core  orbitals,  while  the  QCISD  and  QCISD(T)  used  a  set  of 
frozen  inner  shell  MOs.  All  the  virtual  orbitals  were  included  for  all  the  correlation  calculations. 

In  the  previous  study,  stationary  points  on  the  H  +  OCS  ground-state  PES  were  first  located  at  the 
ROHF  level  (Amos  and  Rice  1992)  using  the  6-31G**  basis  (Hehre,  Ditchfield,  and  Pople  1972; 
Hariharan  and  Pople  1973;  Gordon  1980).  Similar  structures  were  t^en  located  for  each  critical  point  at 
the  correlated  MP2  (M0ller  and  Plessct  1934)  levels  using  UHF  orbttais  (UMP2)  as  the  zeroth-order 
wavefunctions.  One  additional  structure  (structure  [i]  in  Figure  1  of  Rice,  Carlland,  and  Oiabalowski 
[1993])  was  found  at  the  MP2  but  not  at  the  ROHF  level.  Projected  UMP4  energies  were  calculated  for 
each  stationary  point  located  at  the  ROHF  and  UMP2  levels.  UMP2  is  applicable  to  systems  well 
represented  by  tlie  UHF  determinant  (Simandiras  et  al.  1988).  Using  spin  contamination  as  a  metric,  all 
but  one  point  calculated  with  the  6-3 IG**  basis  show  unprojcctcd  UHF  values  of  the  spin  operator  to 
be  0.85  or  less,  the  exception  being  structure  [t]  in  Figure  1  of  Rice,  Cartland,  and  Chabalowski  (1993), 
which  has  =  0.90.  Tliis  suggests  that  the  UHF/6-31G**  wavefunctions  are  reasonable  zeroth-order 
approximations  to  the  wavefunctions. 

In  this  study,  ihc  UMP2  optimizations  were  rerun  for  each  critical  point  located  in  the  previous  study 
using  the  6-311+G(2df,2p)  basis  set  (Krishnan  ct  al.  1980;  McLean  and  Chandler  1980).  Starting 
geometries  for  the  optimizations  were  the  UMP2/5-31G**  structures.  All  the  structures  located  at  the 
UMP2/6-31G**  level  exist  at  the  UMP2/6  31 1+G(2df,2p)  level  and  are  shown  in  Figure  1.  In  addition, 
a  transition  state  connecting  the  di-HSCO  and  cis-HOCS  minima  was  located  (structure  [v].  Figure  1). 
Harmonic  vibrational  frequencies  were  calculated  for  each  stationary  point,  providing  tire  zero-point 
energies  and  characterizing  each  extremum.  UHF  spin  operators  and  zero-point  energies  arc  shown  in 
Table  1.  As  in  the  previous  study  (Rice,  Carilaiid,  and  Chabalowski  1993),  die  stnicture  with  the  largest 
value  of  0.89  corresponds  to  structure  [l]  in  Figure  1.  Tlie  UMP2/6-311+G(2df,2p)  structures  were 
treated  witii  higher  level  correlation  techniques  including  full  fourth-order  (SDTQ)  PUMP4  (Krishnan  and 
Pople  1978;  Krishnan,  Frisch,  and  Pople  1980)  and  QCISD(T)  (Pople,  Head-Gordon,  and  Raghavachari 
1987)  to  generate  single-point  energies  for  all  points  on  the  PES.  'Phese  PUMP4  and  QCISD(T)  results 
will  be  called  "refined  energies."  The  QCISE>(T)AlMP2/6-31  l-t-G(2df,2p)  results  arc  probably  the  highest 
quality  calculations  foi  which  all  the  points  on  the  PES  have  been  studied  and  will  be  called  the  "reference 
set."  Absolute  energies  for  the  optimized  structures  arc  given  in  Table  2.  The  energies  of  these  structures 
relative  to  die  reactants  arc  lisicd  in  Tabic  3  and  illustrated  in  Figure  2. 


i 

I 

1. 1 

[i 


3 


J 


(a) 

(b) 

(c) 

(d) 

(e) 


1.1648  IJ586 

o - c - s 

(1.1536)  (1.5679) 

(1.1642)  (1.5615) 


c-^o 

(1.1296) 

(1.1343) 

S 

(1.5393) 

(15364) 


1.2060 
124.8 


S-L33^H 

(1.3425)** 

(1.3346) 

0.9659 

O - H 

(0.9687) 

(0.9664) 


123.6 


1.7567 


1.1006 

H 


157^  ?9741 

\m.5  iioy 

cVTcdrO 


(i) 


(j) 


1.2964 


1.5303  _  2.1484  0.9735 

-C— - H - O 


(1) 


1.1888 
(1.1845) 

11.1882)  _  ,  ^ 

”^73^^^^^^^^  ^5981 

(173  9)  (1.6125) 

(173.7)  (1.6014) 

1.5748 
(1.5782) 


1.3942 

(1.4007) 

(1.3895) 


(m) 


(n) 


(o) 


15097, 

s' 


1.4167 
Cv.  1.1794 


(P) 


1.7215  V— >sl.‘ 
144.8  O 


1.21'  2799 


O- 


(q) 


(r) 


1.2759  V_7\1.5705 

150.0  ^5 


T(OCSH)=82.0 

»8;y^.3374 
1.8889  ^ 


<S) 


(t) 


1.3336 


1,3375 


(u) 


(v) 


1.2489 


1.6633 


Figure  1.  Optimized  structures  and  geometric  parameters  at  the  UM]*2/6-31  l4-G(2df.2D)  level.  The 
values  of  the  geometric  parameters  of  structures  optimized  at  the  QCISO/6-31 1+G(2df^p) 
and  UMP2/aug-cc-pVlZ  levels  aie  given  in  round  and  square  brackets,  respectively. 
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Table  1.  Zero-Point  Energies  (kcal/mol)  and  of  Species  on  the  H  +  OCS  PES 
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i\s  a  further  check  on  the  basis  set  dependence,  the  reactants,  products,  and  the  two  saddle-point 
structures  labeled  (k)  and  [1]  in  Figure  1  were  reoptimized  at  the  UMP2  level  using  Duitning’s  basis  set 
aug-cc-pVTZ  (Dunning  1989;  Kendall,  Dunning,  and  Harrison  1992;  Woon  and  Dunning  1993).  These 
structures,  and  the  aug-cc-pVTZ  basis,  were  then  used  to  calculate  the  refined  energies  with  full  fourth- 
order  PUMP4  and  QC1SD(T).  Only  these  two  saddle  points  (which  correspond  to  the  two  entrance 
channel  barriers  to  HSCO  formation)  were  studied  due  to  the  extensive  computational  resources  needed 
to  do  these  calculations.  The  aug-cc-pVTZ  and  6-31  l-t-G(2df,2p)  basis  sets  differ  primarily  in  that  the 
former  was  developed  to  optimize  correlation  corrections  to  the  energies  on  atoms  while  the  latter  basis 
was  developed  at  the  11 F  level.  With  respect  to  the  number  and  types  of  basis  functions,  the  sulfur, 
carbon,  atkl  oxygen  atoms  each  have  an  additional  single  primitive  Gaussian  polarization  function  added 
at  both  the  d-  and  f-orbiul  levels  in  the  aug-cc-pVTZ  compared  to  the  6-31 1+G(2df,2p).  In  addition, 
Dunning’s  hydrogen  basis  contains  d-type  polarization  functions  that  do  not  exist  in  the  other  basis 
(Dunning  1989;  Kendall,  Dunning,  and  Harrison  1992;  Woon  and  Durming  1993).  The  valence  spaces 
Ibr  these  two  basis  sets  have  tlie  same  number  of  contacted  AOs  on  hydrogen  and  the  first-row  atoms. 
For  sulfur,  the  6-31 1+Cl(2df,2p)  contains  one  additional  valence  s-  and  p-type  AO,  but  lacks  the  diffuse 
d-  and  l-typc  AOs.  The  aug-cc-pVTZ  contains  165  contracted  AOs  for  H  +  OCS,  while  the 
6-31 1+G(2df,2p)  basis  has  119  AOs,  witJi  Uic  primary  difference  being  attributed  to  the  extra  diffuse 
polarization  functions  in  the  aug-cc-pVTZ  basis.  It  might  be  worth  pointing  out  that  the  6-31 1+G(2df,2p) 
basis  set  is  larger  in  Ixtih  Uie  valence  and  polarization  AO  spaces  than  the  6-3110*  basis  used  in  Poplc 
et  al.’s  G1  and  G2  approaches  (Poplc  ct  al.  1989;  Curti.ss  ct  al.  1991).  Finally,  to  check  what  effect 
conelation  has  on  dte  optimized  structures,  tlic  geometries  for  tiicsc  five  points  were  reoptimized  using 
QCISD/6-31  l-fG(2dl,2i)),  followed  by  an  energy  refinement  at  QCl,SD(T)/6-31 1+G(2df,2p).  The  results 
of  the  extended  treatments  witli  tlic  Dunning  bxsis  set  and  tlie  QCISD  geomcii^  optimizations  arc  listed 
in  Tables  1-4. 

lire  knowledgeable  reader  might  have  noticed  strong  similarities  between  the  current  computational 
approach  and  the  well-known  "Gi"  and  "G2"  approaches  of  Poplc  and  coworkers  (Poplc  ci  al.  1989; 
Curtiss  cl  al.  1991).  We  stress  that  we  have  not  followed  the  Gl  and  G2  guidelines  exactly;  tlicrefore, 
tlicsc  arc  not  G1-G2  type  calculations.  The  basis  sets  used  in  our  calculations  contain  at  most  a  (2d0  set 
of  polarization  functions  lor  nonhydrogcnic  atoms  and  not  the  (3dO  set  recommended  in  tlie  G2  approach. 
This  choice  was  forced  by  computational  limitations.  However,  our  suticiurcs  were  found  at  MP2  (and 
QCISD)  using  a  signincanlly  larger  basis  set  than  recommended  in  Gl  and  G2.  This  was 
done  to  maximize  our  chances  for  accurately  calculating  transition  state  structures  with  elongated  bond 
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arid  allowing  for  jhc  possibility  of  molecular  complexes  to  occur.  In  addition,  no  "higher  level  correction" 
(HLC)  (Poplc  et  al.  1989;  Curtiss  et  al.  1991)  was  applied. 

3.  RESULTS  AND  DISCUSSION 

All  critical  points  found  in  our  previous  study  (Rice,  Cartland,  and  Chabalowski  1993)  were  found 
in  these  larger  calculations.  Additionally,  another  saddle  point  that  connects  the  cis-HOCS  and  c»-HSCO 
wells  was  located.  Tiiis  structure  Oabclcd  (v]  in  Figure  1)  affirms  that  hydrogen  migration  is  possible  for 
this  system,  as  previously  suggested  (BOhmer.  Mikliaylichcnko,  and  Wiuig  1993;  Hduslcr.  Rice,  and  WitUg 
1987).  However,  reaction  through  this  patltway  is  energetically  unfavored,  due  to  tlie  siye  of  the  migration 
barrier.  Isomerization  from  ct,y-HOCS  to  cw-HSCO  at  the  QCISD(T)AJMP2/6-31 1+G(2df,?.p)  level 
required  38.8  kcal/mol.  For  complctcttc.ss,  v'e  sought  and  found  structure  [v]  at  the  UMP2/6-31G**  level 
used  in  our  previous  work.  No  attempt  was  made  to  locate  structure  [v]  at  the  ROHF  level. 

The  structures  obtained  from  Lhe  UMP2/6-31 1+G(2df,2p)  calculations  arc  quite  similar  to  those  found 
with  tile  smaller  6-31G**  basis  in  our  previous  woric  (see  Figure  1  of  Rice,  Cartland,  and  Chabalowski 
[  19931),  showing  a  stability  in  structure  witli  rc.s^icct  to  basis  set  size.  This  convergence  of  suaicturcs  with 
respect  to  basis  set  size  is  furtlier  substantiated  by  comparison  of  the  UMP2/aug-cc-pVTZ  results  given 
in  Figure  1  for  structures  (a],  [bj,  [c],  [k],  and  [1],  represented  by  the  numbers  in  square  brackets. 
Comparison  with  the  nonbracketed  UMP2/6-3 1 1+G(2df,2p)  values  shows  very  little  change  between  the 
bond  lengths  and  angles  calculated  with  tlicse  two  bases.  A  particularly  interesting  fact  is  that  the 
structures  first  located  at  the  ROHF  level  (Rice,  Cartland,  and  Chabalowski  1993)  for  all  tlic  critical  points 
on  the  PES  except  structure  [i]  arc  very  similar  to  those  found  at  the  UMP2  level  for  all  bases  considered 
here. 


We  can  also  compare  the  effect  of  higher  order  correlation  on  structures  [a],  lb],  [c],  [k],  and  [1]  by 

Mt  uiv  Vuiuv'iS  lOi  lii^  i  k  ttJiu  lai  luuiiu  u«aCK.C'U> 

in  Figure  1)  with  the  UMP2/6-31 1+G(2df,2p)  values  (unbrackeied).  The  only  significant  variation  is  the 
S-H  bond  length  in  the  saddle-point  structure  [Ij.  The  UMP2/6-311+G(2df,2p)  and  UMP2./aug-cc-pVTZ 
calculations  predict  this  distance  to  be  1.69  A  and  1.71  A,  respectively.  The  QCISD/6-311+G(2df,2p) 
value,  1.85  A,  is  8-9%  longer  than  the  UMP2  values.  Apart  from  this  S-H  bond  length,  the  structures 
calculated  at  the  ROHF,  UMP2,  and  QCISD  levels  arc  similar,  suggesting  there  will  be  no  gross  structural 
changes  if  higher  levels  of  correlation  arc  used  to  determine  geometries. 
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Tabic  3  gives  a  summary  of  the  energies  for  all  stationary  points  on  the  PES  relative  to  the  H  +  OCS 
asymptote.  The  values  in  Table  3  arc  uncorrected  and  come  directly  from  differences  in  calculated  total 
energies.  The  features  of  die  PES  arc  illustrated  in  Figure  2  and  are  very  similar  to  the  features  of  the 
PUMP4/llMI'2/6-31G**  PES  (see  Figure  2  of  Rice,  Cardand,  and  Chabalowski  [1993]).  The  relative 
barrier  heights  and  well  depths  of  the  PUMP4/UMP2/6-31G**  PES  (Rice,  Cardand,  and  Chabalowski 
1993)  provided  an  cxplanadon  for  the  observed  product  ciicrgy  distributions  for  both  (I)  and  (II). 

For  Reaction  (II),  nonstadstical  product  energy  distributions  were  ob.served.  CO  vibrational  and 
rotadonal  distributions  and  die  SD  rotadonal  disuibuUons  vcrc  colder  dian  predicted  by  statistical  theory 
(BOhmer,  Mikhaylichenko,  and  Wittig  1993;  Nickolaisen  and  Cardand  1993;  Hduslcr,  Rice,  and  Wittig 
1987).  Nickolaisen  and  Cardand  detennined  that  as  much  as  49%  of  available  product  energy  was  in  SH 
vibration  (Nickolaisen  and  Cartland  1993).  The  cxplanadon  for  dicsc  observations  put  fordi  in  our 
previous  study  (Rice,  Caiiland,  and  Cliabalowski  1993)  states  that  upon  fomiation  of  HSCO,  the  molecule 
has  sufficient  energy  to  easily  go  on  to  product,  due  to  the  sizes  of  dte  entrance  and  exit  barriers  and  well 
dcpdi  of  fISCO.  Most  of  diis  excess  energy  will  initially  be  localized  in  the  nascent  S-H  bond  and 
adjacent  I  ISC  bend.  3he  small  exit  ban  ier  and  excess  energy  almost  easurc  that  the  HSCO  molccuie  will 
be  short-lived,  decomposing  before  intramolecular  vibrational  relaxation  (IVR)  is  complete.  Incomplete 
IVR  would  result  in  nonstatistical  product  energy  di.stribudons. 

The  large  exit  channel  barrier  (clo.se  to  die  cndolhcnnicily  of  reaction)  for  Reaction  (I),  on  die  other 
hand,  almost  guarantees  a  long-lived  I  IOCS  complex  before  decomposition  to  product.  The  longer 
lifetime  of  die  molccuie  would  allow  complete  IVR  before  decomposition,  resulting  in  a  statistical 
disuibution  of  product  energies,  as  obseivcd  in  experiment,  lire  qualitative  nature  of  the  surface  has  not 
changed  in  this  study;  therefore,  die  rationalizadon  of  the  product  energy  distributions  for  both  reactions 
offered  in  our  previous  study  (Rice,  Cartland,  and  Chabalowski  1993)  still  hold. 

The  experimental  and  calculated  enthalpies  (corrected  to  T  =  298  K)  for  Reactions  (!)  and  (II)  arc 
given  in  Table  5.  The  dieorctical  numbers  include  zero-point  energy  corrections.  Table  5  contains  results 
for  die  two  basis  sets  6-31 1+G(2df,2p)  and  aug-cc-p VTZ  where  structure  optimizations  were  done  using 
either  the  UMP2  or  QCISD  levels.  All  the  final  energy  refinements  reixirted  in  Table  5  were  done  at  the 
QCISD(T)  level. 
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Tabic  5.  Tcmixiraturc  Corrected  (T  =  298  K)  Theoretical  and  Experimental 
Reaction  Enthalpies  (kcal/ntol) 


Rx.  (I) 

Rx.  (li) 

Experiment: 

AH 

Ref. 

AH 

Ref. 

45.3 

a 

-9.8 

f 

47.1 

b 

-9.8 

gJi 

53.42 

c 

-10.3  ±3.0 

d 

55  ±3.1 

d 

-10.45 

b 

57.2  ±6.0 

-11.0 

i 

-11.3 

a 

-12.1  ±1.2 

c 

■ 

-13 

c 

Theoretical: 

■ 

(1)  AE:  QCISD<T)/6-311+G(2df,2p)/AJMP2/6-311+G(2df,2p) 

59.8 

■ 

-8.8 

Frequencies:  UMP2/6-3 1 1  +G(2df,2p) 

■ 

(2)  AE:  QCISD(T)/6-31l4.G(2df,2p)//QCISD/6-311+G(2df,2p) 

59.7 

■ 

-8.8 

Frequencies:  QCISD/6-311+G(2df2p) 

■ 

(3)  AE:  (X:iSD(T)/aug<c-pVT7y/UMP2/aug-cc-pV'l’Z 

57.3 

■ 

-11.5 

Frequencies:  UMP2/aug-cc-pVTZ 

■ 

*  Nickolaisen  (1991)  ^  Tsunashima  cl  al.  (1975) 

Wagnian  cl  al.  (1982)  *  OldCTShaw  and  Poncr  (1972) 

'  Herzberg  (1945)  Oldcrshaw  and  Smith  (1978) 

“  Hkuslcr,  Rice,  and  Wiilig  (1987)  ‘  Lee,  Siief,  ami  Timmons  (1977) 

*  Chase  ei  al.  (1985) 

3.1  Reaction  (I):  H  +  PCS  OH  +  CS.  Comparison  of  calculated  and  experimental  enthalpies  arc 
presented  in  Table  5.  There  arc  three  sets  of  calculations  listed  under  tlic  theoretical  subheading,  showing 
the  level  of  tlieory  and  basis  set  used. 

Entries  (1)  and  (2)  under  the  Theoretical  subheading  in  Table  5  show  that  the  reaction  enthalpy  for 
(1)  is  the  same  whctlicr  the  geometry  is  optimized  at  the  UMI  2  or  QCISD  levels.  This  simply  verifies 
that  the  characteristics  of  the  products  and  reactant  species  arc  quite  similar  when  determined  by  the  two 
conelation  methods.  Entry  (3)  shows  the  effect  of  the  larger,  "correlation  consistent"  (Dunning  1989-, 
Kendall,  Dunning,  and  Harri.son  1992;  Woon  and  Dunning  1993)  aug-cc-pVTZ  basis  set.  In  this 
calculation,  AH  decreases  by  2.5  kcal/mol  compared  to  the  values  calculated  using  the  6-311+G(2df,2p) 
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basis.  There  is  a  large  range  of  values  given  by  various  experiments  (45.3  to  57.2  kcal/mol)  (Hliusler, 
Rice,  and  Wittig  1987;  Her/.b;rg  1945;  Huber  and  Herzberg  1979;  Nickolaisen  1991;  Wagman  ct  al.  1982; 
Chase  ei  al.  1985)  for  this  reaction.  It  is  difficult  for  us  to  unequivocally  choose  the  "best"  value  from 
this  range,  but  our  calculations  support  the  value  of  Qiase  et  al.  (1985)  (57.2  ±6.0  kcal/mol).  All 
calculated  AH’s  in  Table  5  fall  within  the  experimental  uncertainty  of  the  Chase  et  al.  (1985)  value,  with 
the  QCISD(T)/aug-cc-pVTZ  value  agreeing  most  closely.  Such  close  agreement  for  entry  (3)  is  most 
likely  fortuitous,  but  it  docs  lend  support  to  the  possible  quantitative  nature  of  these  calculations. 

Figure  2  shows  that  the  paths  leading  to  formation  of  the  products  of  Reaction  (1)  [as  well  as  for 
Reaction  (II)]  arc  anytliing  but  simple.  Possible  direct  pathways  to  formation  of  the  products  of  (I)  are: 


[m]  [u] 

H  +  OCS  _ k  trans  -  HOCS  _ ^  OH--CS  _ ^  OH  +  CS 


(I.A) 


[n] 

H  +  OCS  _ K  cis  -  HOCS  _ k  OH  ^  CS 


(I.B) 


Also,  some  reactions  involving  isomerization  of  tlic  four-body  species  that  could  lead  to  products  are: 


Ijl  [p]  [u] 

H  +  OCS  _ kHOCS  _ k  rra/ur-HOCS  _ k  OH  — CS 


OH  f  CS 


(I.C) 


In]  [t]  [u] 

H  +  OCS  _ kcis-KOCS  _ k  rrans-HOCS  _ k  OH--CS 


OH  +  CS 


(I.D) 


[m]  [t] 

H  +  OCS  _ kira/is-HOCS  _ k  cts-HOCS  _ k  OH  +  CS 


(I.E) 


P]  [v] 

H  +  OCS  _ k  cis  -  HSCO  _ k  cis  -  HOCS  _ k  OH  +  CS 


(I.E) 
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PI  [v]  m 

H  +  OCS  _ ^  cis-HSCO _ ^  cis-HOCS  _ ^  iriUu-HOCS 


[u] 


OH  — CS 


OH  +  CS 


(I.G) 


We  did  not  include  all  isomerization  reactions  that  involve  fonnation  of  c/s-HSCO  (through  formation  of 
trans-USCO  either  directly  or  through  foimation  of  HCOS):  the  few  listed  above  indicate  the  possible 
pathways  that  can  be  accessed  given  enough  energy. 

Although  we  cannot  conclude  that  the  most  probable  reaction  paths  are  (I.  A)  and  (I.B)  without  doing 
dynamics,  they  arc  the  most  direct.  The  two  entrance  channel  barriers,  [m]  and  [n],  for  these  direct  paths 
arc  rubstaniially  lower  (<25  kcal/mol)  than  the  reaction  cndothermicity  (>50  kcalAnol).  Tlie  structures 
suggest  that  hydrogen  prefers  a  sidc-on  attack,  with  C-O-H  angles  of  1 1 8“  and  1 12°  for  the  cis-  and  trans- 
spccies,  respectively.  Wc  were  unable  to  locate,  as  in  the  previous  study  (Rice,  Cartland,  and  Chabalowski 
1993),  an  exit  channel  barrier  leading  from  c«-HOCS  to  product.  However,  the  path  leading  from  trans- 
HOCS  must  overcome  a  barrier,  [u],  which  lies  below,  but  close  to  the  endothcrmicity  of  reaction.  To 
pass  from  [u]  to  products,  the  reaction  goes  through  an  interesting  intermediate  that  is  stable  (with  respect 
to  [u])  by  only  0.4  kcal/mol.  This  stable  intermediate,  [i],  represents  a  linear  OH--CS  geometry  when 
hydrogen  has  now  rotated  to  point  toward  carbon  at  a  large  separation  of  2.15  A. 

Tire  barrier  to  formation  of  tlic  HCOS  species  (hydrogen  is  attached  to  the  central  carbon  atom)  is  less 
than  half  of  tliosc  leading  to  direct  HOCS  fomiation.  However,  to  go  on  to  tfie  product  of  (I),  HCOS  must 
overcome  a  large  isomerization  barrier,  [p],  to  form  rra^w-HOCS.  This  isomerization  requires 
approximately  16  kcal/mol  more  to  reach  rra/u-HOCS  [g]  than  docs  the  direct  attack  of  H  on  O.  Also, 
tire  barrier  to  formation  of  ds-HSCO  is  the  lowest  barrier  on  this  PES;  however,  the  isomerization  to  cis- 
HOCS,  necessary  to  go  to  the  products  of  (1),  is  substantial  (~50  kcalAnol).  Since  the  energy  of  the 
products  is  greater  than  all  barriers  on  tliis  PES,  all  of  the  paths  listed  above  arc  energetically  accessible 
for  reaction  conditions  that  will  lead  to  (I). 

3.2  Reaction  (ID:  H  +  OCS  SH  +  CO.  The  experimentally  determined  reaction  enthalpies  for 
(II)  range  from  -9.8  to  -13  kcal/mol  (Hduslcr,  Rice,  and  Wittig  1987;  T.sunashima  et  al.  1975;  Lee,  Stief, 
and  Timmons  1977;  Herzberg  1945;  Huber  and  Herzberg  1979;  Nickolai.scn  1991;  Wagman  et  al.  1982; 
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Chase  et  al.  1985;  Oldcrsi.aw  and  Porter  1972;  Oldershaw  and  Smith  1978).  If  one  assumes,  as  we  did 
in  our  comparisons  for  Reaction  (I),  that  the  critically  evaluated  enthalpy  given  by  Chase  et  al.  (1985)  is 
the  most  reliable  number,  then  experimental  reaction  enthalpy  for  Reaction  (II)  is  -12.1  ±1.2  kcal/mol. 
Table  5  shows  that  UMP2  and  QCISD  reaction  enthalpies  calculated  with  the  6-31 1+G(2df,2p)  basis  arc 
predicted  to  be  -8.8  kcal/mol,  or  3.3  kcal/mol  higher  than  the  Chase  ct  al.  (1985)  value  and  1.0  kcal/mol 
above  the  highest  reported  experimental  value  (Tsunashima  et  al.  1975;  Oldershaw  and  Porter  1972; 
Oldershaw  and  Smith  1978).  However,  the  results  from  the  larger  "correlation  consistent"  basis  set 
calculation  [entry  (3),  Table  5]  give  AH  =  -11.5  kcaljtaol,  well  within  the  ±1.2  kcal/mcl  uncertainty  in 
the  Chase  et  al.  (1985)  value.  The  results  of  AH  calculated  for  Reactioris  (I)  and  (II)  show  that  the 
energetics  of  tltis  system  are  not  converged  with  respect  to  basis  set  at  the  6-311+G(2df,2p)  level,  but 
approach  basis  set  convergence  with  the  aug-cc-pVTZ  basis,  at  least  if  one  defines  convergence  as 
predicting  reaction  enthalpies  that  fall  witliin  experimental  uncertainty.  It  must  be  pointed  out  that  moving 
from  the  6-31 1+G(2df,2p)  to  tlic  aug-cc-pVTZ  basis  is  accompanied  with  a  sizable  increase  in 
computational  effort.  Given  the  range  of  experimental  data  for  (II),  and  that  the  Chase  et  al.  (1985)  value 
is  at  the  low  end  of  this  range,  it  is  probably  safe  to  say  that  the  3.3  kcal/mol  eiTor  obtained  using  the 
6-311+G(2df,2p)  basis  is  an  upper  limit,  and  that  the  actual  caor  is  less. 

The  portion  of  the  PES  for  Reaction  (II)  is,  as  for  Reaction  (I),  very  complex.  Figure  2  shows  several 
pathways  for  reaching  the  products  SH  +  CO,  many  involving  isomcrizations.  Some  of  them  include: 


M  [r] 

H  +  CX:S  _ ^  trans  -  HSCO  _ ^  SH  +  CO 


(II.A) 


D]  [s] 

H  +  OCS  _ ^cis-HSCO  _ ^SH  +  CO 


(II.B) 


Ul  lo]  W 

H  +  OCS  _ w  HCOS  _ ^  rrowf-HSCO  _ v  SH  +  CO 


(II.C) 


m  [ql  [r] 

H  +  OCS  _ V.  cis-HSCO  _ ^  rrans-HSCO  _ ^  SH  +  CO 


(II.D) 


[k]  [q]  [s] 

H  +  OCS  _ ^  trans  -  HSCO  _ w  Ci.^  -  HSCO  _ ^  SH  +  CO 


(Il.E) 
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SH  +  CO 


[n]  [v]  [s] 

H  +  OCS  _^cw-HSCO  cis-HSCO  y 


(II.F) 


As  in  our  analysis  of  paths  resulting  in  Reaction  (I),  we  have  not  included  all  of  the  possible  isomerization 
reactions  (such  as  formation  of  trans-HOCS,  then  isomerization  to  cij-HOCS,  etc.)  in  the  previously 
mentioned  list.  However,  they  are  possible  under  the  correct  experimental  conditions. 

Reactioas  (II.Q  and  (II.F)  are  high-energy  paths,  and  probably  minor  contributoa  to  the  reaction. 
High-energy  pathway  (li.F),  and  derivatives  of  it.  appeared  with  the  discovery  of  the  c«-HOCS  to  cis- 
HSCO  hydrogen-migration  barrier  [v].  This  transition  state  is  the  highest  energy  barrier  along  any  of  the 
pathways  involving  H-atom  migration  (~40  kcal/mol).  Once  over  this  barrier,  the  system  has  ample 
energy  to  proceed  to  SH  4-  CO.  Any  paths  leading  to  cis-HOCS  followed  by  H-atom  migration  to  the 
sulfur  could  represent  the  hydrogen-migration  patii  suggested  previously  as  a  possible  mechanism 
(Btihmer,  Mikhaylichenko,  and  Wittig  1993;  HSusler,  Rice,  and  Wittig  1987).  In  the  more  recent  of  these 
two  experimental  papcis  (Bdhmer,  Mikhaylichenko,  and  Wittig  1993),  the  authors  used  sUuctuies  and 
relative  energies  from  our  earlier  study  (Rice,  Cartland,  and  Chabalowski  1993)  to  cast  doubt  on  this  path 
as  a  viable  mechanism  for  producing  the  observed  vibrationally  coo!  CO  product  They  also  argued  that 
a  sudden  shortening  of  the  elongated  C-0  bond  in  HOCS  should  occur  (e.g.,  structure  [h])  when  H 
undergoes  a  1.3  shift  from  O  to  S  and  or.  to  products.  This  should  produce  a  vibrationally  excited  CO 
stretch  thus  giving  an  energy  distribution  considerably  hotter  in  the  stretch  mode,  just  the  opposite  to  what 
is  observed  experimentally  (BOhmer,  Mikhaylichenko,  and  Wittig  1993:  Hausler,  Rice,  and  Wittig  1987). 
Tlie  current  results  further  support  their  interpretations. 

High-energy  pathway  (II. C)  begins  with  H  attacking  the  central  carbon  atom  to  form  the  stable  HCOS 
intermediate  [d],  discussed  in  the  previous  section.  From  here  the  H-atom  migrates  to  sulfur  over  the 
largest  barrier  along  this  path,  [o].  Isomerization  from  HCOS  to  rranj-HSCO  requires  ~30  kcal/mol. 
i*fter  crossing  [o],  the  nascent  rra«i-HSCO  will  have  sufficient  excess  energy  to  easily  go  onto  products. 

One  of  the  more  interesting  results  from  this  work  comes  from  calculations  of  the  barriers  for  the  two 
low-energy  pathways  for  (II)  at  higher  levels  of  theory.  The  pathways  involve  direct  attack  of  the  H  atom 
on  the  sulfur.  In  the  previous  study  (Rice,  Cartland,  and  Chabalowski  1S>93),  the  lowest  entrance  channel 
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barrier  for  this  reaction  was  formation  of  ds-HSCO  (12  kcalAtxol);  the  next  higher  barrier  (formation  of 
rrani-HSCO)  was  3  kcal/mol  higher.  This  energy  separation  between  the  two  barriers  showed  that  the 
major  contribution  to  the  thermal  rate  would  be  due  to  reaction  via  the  ds-HSCO  formation.  In  this  study, 
however,  the  barrier  to  formation  of  rra/u-HSCO  is  more  strongly  influenced  with  the  increase  in 
correlation  and  basis  set.  The  barriers  {[k]  (rro/w-HSCO)  and  [[1]  (ds-HSCO))  have  calculated  heights 
of  7.1  kcal/mol  and  6.3  kcal/mol  at  the  QCISD(T)//UMP2/6-31 1+G(2df,2p)  level,  respectively.  As  seen 
in  Figure  2,  once  the  entrance  barriers  [k]  and  [1]  are  overcome,  the  molecules  have  enough  excess  energy 
to  overcome  the  exit  barriers  [r]  and  [s]  {trans-  and  cis-,  respectively)  and  go  on  to  products,  making  tlie 
entrance  channel  barriers  the  rate  limiting  step  along  each  reaction  path. 

Having  alluded  tc  the  probable  uncertainties  (supra  vide)  in  the  present  calculations,  the  energy 
difference  of  0.8  kcal/mol  between  tlicsc  two  barriers  is  too  small  to  make  any  strong  statement  about  the 
relative  ordering  of  the  two  barriers.  We  carried  out  a  scries  of  calculations  for  these  two  entrance  channel 
barriers,  in  an  attempt  to  quantify  our  results.  The  barrier  heights  for  [k]  and  [1],  respectively,  arc 
7.0  kcal/mol  and  6.8  kcal/mol  at  the  QCISD(T)//QCISD/6-311+G(2df,2p)  level,  and  5.3  kcal/mol  and 
5.2  kcal/mol  at  the  QCISD(T)//UMP2/aug-cc-pVT':  level.  From  what  we  consider  our  best  calculations 
(QCISD(T)//UMP2/aug-cc-pVTZ),  the  entrance  channel  barrier  to  formation  of  cw-HSCO  lies  only 
0.10  kcal/mol  below  tlic  barrier  to  formation  of  trans-HSCO. 

One  notable  exception  to  energetic  ordering  of  barriers  to  formation  of  cis-  and  fr,’.»a-HSCO  occurs 
in  tlic  QCISD  calculations.  The  QCISD  without  triples  reverses  the  order  with  the  (rans-  barrier  now 
lying  ~l-2  kcal/mol  below  iJic  cis-  barrier.  The  inclusion  of  triples  puts  the  cis-  barrier  lower  than  the 
trans-  barrier,  as  seen  in  all  of  Uic  oilier  PUMP4  and  QCISD(T)  results.  Apparently,  the  triples  play  an 
important  role  in  the  details  of  the  relative  energetics  of  these  transition  states. 

It  is  also  interesting  to  compare  the  PUMP4/aug-cc-pVTZ  barrier  height  witli  the 
QCISD(T)/aug-cc-pVTZ  value,  llic  PUMP4  predicts  a  1.9  kcal/mol  difference  between  the  two  barriers, 
notably  larger  than  tlie  0.10  kcal/mol  difference  predicted  by  QCISD(T).  We  sec  the  same  trend  ol 
QCISD(T)  narrowing  the  energy  difference  between  the  cis-  and  trans-  barriers  compared  to  the  PUMP4 
separation  for  the  6-31 1+G(2df,2p)  basis.  This  could  attribute  to  the  "often  unsatisfactory"  convergence 
of  the  MP  scries  (Kri.shnan,  Frisch,  and  Poplc  1980). 
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Figure  3  is  an  attempt  to  detect  a  trend  in  relative  energies  for  [k]  and  [1]  with  respect  to  the  level  of 
correlation  correction  and  basis  sets  used  in  the  calculations.  Admittedly,  we  cannot  state  with  certainty 
that  the  results  in  any  one  column  of  Table  3  arc  more  accurate  than  another.  However,  based  on  our 
conclusions  regarding  reaction  enthalpies,  it  will  be  assumed  that  the  calculations  become  more  accurate 
as  one  proceeds  from  left  to  right  in  Figure  3  and  in  Table  3.  An  approach  to  convergence  for  the  barrier 
heights  of  [k]  and  [1]  with  respect  to  level  of  basis  set  and  correlation  theory  is  apparent  in  Figure  3.  The 
converged  value  for  each  is  in  the  range  of  5-6  kcalAnol. 

Tsunashima  ct  al.  (1975)  report  a  pre-exponential  factor  for  Reaction  (11)  that  is  smaller  than  many 
hydrogen-abstraction  reactions  and  conclude  that  (T)  has  a  low  entropy  of  activation.  They  make  a  rough 
estimate  for  the  entropy  of  the  activated  complex  and  state  that  agreement  witli  experiment  can  be 
achieved  if  the  complex  has  a  degenerate  bending  mode  of  dOO  cm'^  (They  have  assumed  a  linear 
stiucture  for  [k]  or  [1]  and  hence  the  statement  of  a  "degenerate"  bending  mode.)  They  then  suggest  that 
this  "...  indicates  a  rather  tight  stiucturc  for  the  activated  complex."  All  three  levels  of  theory  used  in  this 
work,  as  well  as  in  our  previous  study,  predict  two  nonlinear  transition  states  leading  to  the  cis-  and  trans- 
isomers  of  HSCO.  The  geometric  parameters  shown  in  Figure  1  corresponding  to  [k]  and  [1]  show  the 
OCS  angles  as  173°  and  17 1°  for  the  trans-  and  cis-  barriers,  respectively.  The  HSC  angles  arc  1 03°  and 
115°  for  [k|  and  [1],  respectively.  The  lowest  energy  bending  modes  (from  the  reference  set)  have 
hequencics  of  479  cm"'  and  415  cm'^  for  [k]  and  [1],  mspcctively  (Table  4).  Tltc  UMP2/aug-cc-pVTZ 
results  give  466  cm'^  and  396  cm  ',  respectively  (Table  4).  The  sizes  of  the  vibrational  frequencies  of 
the  bound  modes  of  tlic  saddle  points  suppor*.  tlic  analysis  provided  by  Tsunashima  el  al.  (1975)  of  a  tight 
activated  complex,  aliliough  we  do  not  show  linear  activated  complexes  that  lead  to  products. 

Before  we  can  compare  the  theoretical  energies  to  the  experimental  rate  measurements,  corrections 
for  zero-point  energy  must  be  applied.  Tlicre  ate  no  experimental  frequencies  available  for  transition 
states  [k]  and  [1]  to  use  to  tnakc  these  corrections.  In  addition  to  the  zero-point  energy  correction, 
tunneling  must  be  considered.  Schalz  (1987),  in  his  review  of  the  effect  of  tunneling  on  bimolecular 
collisions,  urges  the  reader  to  use  "extreme  caution"  when  comparing  barriers  (such  as  those  presented 
here)  to  experimentally  determined  activation  energies,  due  to  possibly  significant  tunneling  effects.  He 
lists  several  studies  that  show  substantial  contribution  to  the  rate  due  to  tunneling  for  hydrogen  reactions. 
He  cites  specific  examples  for  reduction  of  activation  energy  due  to  tunneling;  the  reactions  cited  arc 
K  t-  Hj  and  H  +  HD.  Arrhenius  fits  to  the  rates  between  200  K  and  300  K  show  Uiat  the  activation 
energy  is  lowered  by  3.8  kcal/mol  and  3.7  kcal/mol  for  the  H  +  Hj  atid  H  -r  HD  reactions,  respectively. 
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when  tunneling  is  included.  Additionally,  Schatz  also  shows  that  the  pre-exponential  factors  for  these 
reactions  change  by  orders  of  magnitude  when  tunneling  is  included.  It  is  reasonable  to  expect  that 
tunneling  will  contribute  to  the  thermal  rate  for  this  reaction.  Based  on  tlie  work  reviewed  by  Schatz 
(1987),  this  contribution  could  be  as  large  as  3-4  kcal/mol.  Any  further  attempt  to  quantify  this 
contribution,  however,  is  outside  the  scope  of  this  study. 

Kudla,  Schatz.  and  Wagner  (1991),  in  their  study  of  the  OH  +  CO  (the  isovalent  reverse  of  reactions 
presented  here)  considered  the  deray  rate  of  tlic  four-body  HOCO  reaction  intemiediate  to  H  +  CO2.  They 
found,  in  a  tiinneling-corrected  RRKM  calculation,  that  the  rate  of  decay  is  very  sensitive  to  character  of 
the  saddle  point  between  HOCO  and  H  +  CO2.  For  example,  they  showed  that  tunneling  has  little  effect 
on  the  rate  (~33%  increase)  if  the  barrier  has  an  imaginary  frequency  of  5231  cm■^  but  a  large  effect 
(742%  increase)  when  the  frequency  is  1,5001  cm'*.  Tlie  imaginary  frequencies  (at  the 
UMP2/aug-cc-pVTZ  level)  for  tlic  saddle  points  leading  to  as-  and  tra/is-HSCO  are  1,3541  cm'*  and 
123i  cm'*,  respectively.  This  suggests  that  tunneling  might  be  significant  for  cls-HSCO  formation. 

The  zero-point  energy  corrected  values  for  the  QCISD(T)/AJMP2/aug-cc-pVTZ  barriers  to  formation 
of  dv-HSCO  and  iraru-HSCO  are  5.7  kcal/mol  and  8.4  kcal/mol,  respectively.  The  cis-  barrier,  after  zero- 
point  energy  correction,  is  within  1.9  kcal/mol  of  the  reported  experimental  activation  energies,  without 
considering  tunneling.  The  zero-point  energy  corrected  value  for  tlie  trans-  species  indicates  that  in  order 
to  be  competitive  witli  the  cis-  species,  tunneling  must  lower  tire  activation  energy  by  ~4.5  kcal/mol.  If 
this  system  behaves  as  lire  analogous  H  +  COj  reaction  cited  by  Kudla,  Schatz,  and  Wagner  (1991), 
tunneling  would  be  less  for  fonnation  of  rra«.v-HSCO  than  cls-HSCO  due  to  Uic  characterization  of  the 
saddle  point  provided  by  the  imaginary  frequencies.  This  leads  to  the  conclusion  that  reaction  via  the 
formation  of  cli-HSCO  is  probably  tlic  dominant  path  for  thennal  reaction.  It  also  suggests  that  the 
theoretical  results  presented  here  can  be  considered  to  be  in  good  agreement  with  the  experimental  rate 
measurements.  Even  the  presumably  less  accurate  reference  set  energy  for  the  cis-  entrance  barrier 
(including  zero-point  energy  correction)  still  falls  within  3  kcal/mol  of  the  experimental  value,  without 
considering  tunneling. 

Finally,  we  attempt  to  asse.ss  the  uncertainty  in  tire  theoretical  results.  To  do  this  successfully,  one 
must  have  reliable  experimental  numbers.  We  have  chosen  (assumed)  the  critically  evaluated  reaction 
enthalpies  of  Chase  et  al.  (1985),  and  tlic  experimental  activation  energy  of  3.9  kcal/mol  (Tsuna.shima  cl 
al.  1975;  Lee,  Stief,  and  Timmons  1977)  for  Reaction  (11)  to  meet  tliis  requirement.  Tire  reference  set  of 
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calculatiotis  gives  a  reaction  enthalpy  differing  from  experiment  by  2.6  kcaVmol  for  Reaction  (I)  and 
3.3  kcal/inol  for  Reaction  (II).  The  cis-  entrance  barrier  (rcfercnce  set  energy,  witli  zero-point  correction, 
no  tunneling)  is  within  3  kcal/mol  of  the  experimental  activation  energy.  Putting  error  bars  on  every  point 
on  the  PES  is  not  possible,  but  in  light  of  the  agreement  between  theory  and  what  little  we  know 
experimentally  about  [k]  and  [1],  the  ±3.3  kcal/mol  uncertainty  seems  a  reasonable  estimate  for  most  points 
on  the  PES. 

4.  CONCLUSIONS 

We  have  attempted  to  refine  features  of  the  H  +  OCS  potential  energy  surface  from  our  previous 
study,  using  better  basis  sets  and  potentially  better  correlation  techniques.  In  the  previous  work  (Rice, 
Cartland,  and  Chabalowski  1993),  we  located  structures  at  Uic  UMI^/b-SlG*-*  level,  followed  by  energy 
I'cfijicmcnts  using  PUMP4.  All  structures  located  in  tlic  cailicr  work  were  found  in  this  study  at  the 
UMP2/6-31 1+G(2df,2p)  level.  We  found  an  additional  transition  state  s’ructure  corresponding  to  hydrogen 
migration  from  c«-HOCS  to  cL-HSCO. 

The  general  features  of  tlic  PES  (structures,  relative  baniers,  and  local  minima)  arc  qualitatively  the 
same  as  those  in  our  previous  study  (Rice,  Cartland,  and  Chabalowski  1993).  These  features  were  used 
to  explain  the  observed  product  energy  distributions  for  Reactions  (1)  and  (II).  The  nonstatistical  product 
energy  distributions  for  SH  and  CO  could  be  explained  by  the  relative  sizes  of  the  entrance  and  exit 
channel  barriers  and  well  depths  for  tlic  HSCO  molecule.  Upon  crossing  the  larger  entrance  channel 
barrier  into  a  region  of  die  PES  corresponding  to  the  HSCO  minimum,  the  molecule  has  excess  energy 
to  overcome  the  small  exit  channel  banicr,  which  could  easily  occur  before  IVR  is  complete.  Tlie 
statistical  product  energy  distributions  for  OH  +  CS,  on  tlie  other  hand,  can  be  explained  in  light  of  the 
large  exit  channel  barrier  leading  to  OH  +  CS.  The  size  of  the  barrier  is  close  to  the  reaction 
cndothermicity,  which  must  be  localized  in  tlic  reaction  coordinate  before  reaction  can  proceed.  This 
ahnosi  ensures  Uiai  tlie  moiecuie  will  be  long-lived,  resulting  in  more  complete  iVR  and  statistical  product 
energy  distributions.  These  conclusions  remain  unchanged,  as  the  relative  features  of  the  PES  are 
unchanged:  The  larger  barriers  remain  large,  and  the  smaller  barriers  arc  small.  TTic  well  depths  do  not 
decrease  appreciably.  Tire  most  important  differences  in  the  PES  arc  the  decreases  in  the  entrance  channel 
baiiicrs  leading  to  cis-  and  //-aru-HSCO.  The  best  uncorrceted  estimates  from  this  study  predict  tliat  the 
trans-  and  cis-  barriers  are  almost  equal  in  magnitude.  The  best  calculated  theoretical  values  for  the  trans- 
and  cis-  barriers,  excluding  any  quantum  or  dynamic  corrections,  are  5.3  kcal/mol  and  5.2  kcal/molc. 
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respectively.  Upon  applying  zero-point  energy  corrections,  the  energy  separation  between  the  two  barriers 
is  2.65  kcal/mol,  with  the  cis-  barrier  being  only  1.9  kcal/ntol  higher  in  energy  than  the  measured 
activation  energy  (~3.9  kcaVmol).  Making  a  reasonable  approximation  for  tunneling  and  inclusion  of  the 
zero-point  energy  brings  the  value  for  tlic  banter  to  formation  of  c/j-HSCO  into  good  agreement  with  the 
experimental  measurements.  Tlie  zero-point  energy  corrected  barrier  to  formation  of  frazu-HSCO,  on  the 
other  hand,  is  ~4.5  kcal/mol  larger  titan  the  experimentally  measured  activation  energy,  suggesting  that 
unless  tunneling  plays  a  significantly  greater  role  in  overcoming  the  trans-  entrance  barrier  titan  it  docs 
for  the  cis-  barrier,  this  pathway  will  not  be  competitive  with  the  cis-  channel  for  ilicrmal  reaction. 

Additionally,  we  have  located  a  previously  unrcported  saddle  point  between  tlic  cis-HOCS  and  cis- 
H.SCO  species,  which  confimis  the  hydrogen-migration  route  hypothesized  by  Hduslcr,  Rice,  and  Wittig 
(1987).  However,  the  magnitude  of  the  calculated  barrier  suggests  tltis  is  not  an  energetically  favored  path 
to  formation  of  SH  +  CO. 

In  addition  to  addressing  tlic  tliennal  rate  exixirimcnts,  we  have  compared  effects  on  tlic  entrance 
channel  barrier  heights  and  reaction  cntlialpics  from  different  correlation  treatments  and  basis  sets.  'ITic 
calculated  reaction  cntlialpics  tliat  agree  most  closely  with  cxjicrimcnt  are  tlic  QC1SD{T)  calculations  using 
die  aug-cc-pVTZ  basis.  These  same  calculations  provide  the  lowest  entrance  channel  barriers  to  reaction. 


Intentionally  lelt  blank. 
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